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ABSTRACT 



Interaction between an electron beam and a fully ionized plasma has 
been studied with a view towards its application in a structure-less 
traveling wave tube. Three basic approaches, of varying degree of rigor, 
to the problem have been pursued and analytical solutions for the circularly 
symmetric case obtained. Comparisons between the methods of analysis 
are made. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Professor Glenn A. Gray of the U. S. Naval 
Postgraduate School in this investigation. 
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1 . Introduction 



Investigations into wave propagation in plasmas have been conducted 
by a number of workers (1), (2), (3) with a view towards employment in 
plasma diagnostics and traveling wave tubes. Extensive bibliographies are 
contained in these works and the reader is referred to these for previous 
work in the field. 

Basically, the problem is a boundary value problem in cylindrical co- 
ordinates (r,9,z) as depicted in Fig. 1 with a static magnetic field parallel 
to the z-axis. Throughout this paper, the case of coupling to a helix will 
be considered as this is a common method of coupling into or out of a 
traveling wave tube. The mathematical model used for the helix is the 
simplest of those that have been developed, that of the "sheath helix". 

The reader is referred to any standard work, such as Pierce (4) for details 
on this point. The results will be used without comment in this paper. 

The mathematical model used for the plasma will be that customarily 
used in studies of this nature unless specifically stated otherwise. This is 
the "cold" plasma wherein effects of collisions, recombinations, neutrals, 
and thermal motion are ignored and the electrons are considered as forming 
a "cloud" against a background of positive ions which provide overall (DC) 
neutrality of charge, but do not otherwise appreciably contribute to the 
problem due to the relatively large mass of the ions which renders them 
practically stationary. 

Three approaches to the problem are presented and compared. These 
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0 < r < a Electron Beam 
• a < r < b Plasma 

b < r < c Free Space 

i 

r = c Helix 

r>c Free Space 

Figure 1 

Geometry of the Problem 
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are presented in order of increasing complexity, if not rigor, and the results 
are compared using the third method as the standard of comparison 

The first method is an extension of the works of Boyd, Gould and 
Trivelpiece and has, as the major distinguishing features, the assumptions 
that the magnitude of the longitudinal propagation constant is much greater 
than that of free space and that the electric field may be represented as 
derived from a scalar potential. This is Trivelpiece's ’'slow wave" or 
"quasi-static" approximation. This method has the advantage of simplicity 
albeit at the price of rigor, but the extreme simplicity alone is of considerable 
practical value provided, of course, that the results provide a reasonable 
approximation to the true case. A feature of this "slow wave" approximation 
is that a TE mode of propagation is denied by the first assumptions . While 
of no great consequence when dealing with drift tubes or waveguides, this 
is troublesome when one attempts to derive expressions for interactions wiih 
a helix. Matching boundary conditions at the "sheath helix" requires both 
TE and TM types of solutions as is stated in Hutter (5) and can quickly 
be demonstrated. Thus, a dilemma presents itself* In attempting to extend 
this simple method, a free space TE solution will be assumed within the 
electron beam and plasma regions and the inconsistency ignored. 

The second method is an extension of the works of Rigrod and Lewis (6 , 
and Brewer (7), the latter being essentially a generalization of the former. 

This method, as developed for electron beam studies, solves Maxwell's 
equations in a region containing charge, and, through a perturbational 
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approach, takes the effect of the charge into consideration in the boundary- 
value problem by replacing the rippled beam by an equivalent smooth beam 
with a surface current density. Brewer's model for the electron beam is 
simply a beam of electrons and does not postulate positive ion neutrali- 
zation of the beam as some other common analyses do. The plasma, in 
this type of analysis, is treated as the limiting case of an electron beam 
in a plasma with zero charge density. Brewer, in his paper on the subject 
(7), obtains only a TM solution although a TE solution is not negated as 
in Trivelpiece's analysis. In this paper, an approximate TE solution is 
obtained, and interestingly, is shown to be coupled to the TM solution 
in such a manner that, if the TE solution is identically zero, then the TM 
solution is also zero, with the converse not necessarily true. 

The third method of analysis is considered the most rigorous and is the 
most complex, mathematically. It consists of solving Maxwell's equations 
in an anisotropic media using Kales' (8) method of solution. This method 
provides an exact solution of the mathematical model as described above 
and further amplified in appendix A. Interesting features of this method are 
the requirement of coupled TM and TE modes of propagation which cannot 
be zero independently, and a mode degeneracy with non-orthogonal modes. 

In this paper, the subscript o will be used to denote dc (static) 
quantities and the subscript 1, time varying quantities. Unless noted to 
the contrary, all quantities will be assumed to vary as 
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F ~ ft (^> 0 ) i)£ 



ajt. 



(1 1 ) 



with a z-dependence such that 



A£=-ff(^) n.2) 

1u.'t 

Following the usual procedure, the factor e will be understood and not 
written explicitly . 
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2. Method I, Trivelpiece's "Slow Wave" Approximation 



Following Trivelpiece (1), it is assumed that 

VX H/ ~~ ^ 0 (2.1) 

which then allows us to represent the electric field vector as derived from 
a scalar potential, i.e. 

B, - ~ V (2.2) 



From Maxwell's equations, we have 

V-D, = V-a-E,)= v-l -V^ = 0 



where ~ is the tensor dielectric constant given by 
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A simple derivation of equation (2.4) is given as Appendix A although it is 
worth noting, in passing, that the same result for the case of a plasma 
alone may be obtained from the Boltzmann transport equation (10) with far 
fewer restrictions upon the derivation. 

Proceeding formally, one then obtains 



- j- r ^ 

st d/vL 




f 











(2.5) 
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If a product type of solution is assumed, 



(J), — Rfa) ^ {2 ° 6) 

the following differential equation is then obtained 
iA r 0 (2«7) 

the solution of which is 





C ccxifa 9) -hOytU^CnQ) 




( 2 . 8 ) 



where A, B, C, and Dare arbitrary constants. Following the usual pro- 
cedure in waveguide type propagation problems this will be expressed as 



where 




( 2 .-.) 



9. 

T - 



& . 

e,i 



(2.10) 



We then have 

e u = - + e r< (t~) 



i 



( 2 . 11 ) 




( 2 . 12 ) 



7 



{2 A3) 




where the prime indicates differentiation with respect to the total argument 
While it is implicit in the assumption behind equation (2.2) that the 
time varying magnetic field is essentially zero with respect to the electric 
field, the curl H equation 

V K H, — ^ ^ t.* E/ (2.i4) 



will be used to obtain approximate values of the magnetic field components 

H will be set equal to zero, otherwise the existence of another mode of 
z 

propagation would be allowed which would contradict the basic assumption 
equation (2.2). Taking the components of equation (2.14), one obtains 









'// W-'O 



(2 o 15) 




~jr uJe o [ 



(2.16) 




We then obtain (where the primes again denote differentiation with respect 
to the total argument) 
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(2 ] 8 ) 



B h = - AT ri(T~) +8 TAfcpK)) + 

% 4g ^(U) +B»JTs^ J 




■h 



( ^ fA JJT.) +BN„(t4) i| (2 . 19) 

From this point on, only the axially symmetric case (n=0) will be considered. 
Also, only the time varying field quantities are involved so the subscript 1 
will be omitted to simplify the notation. 

£21 - AYJ^(TsC) + £M£( 7>0 (2.20) 



f/L - AVJ-'fa) +■ BT^ (TrC) 




4>jJ£q 



-// 



/a TJ, (7XJ fBTty fuj 



( 2 . 22 ) 
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(* +■£ tn, (t^ 



12 . 23 



Equations (20) through (23) are the field quantities derived from equa’ions 1 
and (14), to be used in the boundary matching problem. 

As was mentioned in the introduction, it is impossible to match a TM 
solution alone to the "sheath helix" using the standard boundary condi- 
tions as given by Pierce (4) or Beck (9). Some form of the TE solution must 
also be used. For this analysis, a free space TE solution will be assumed 
with more to be said upon this assumption later It is to be noted that up 
to this point, no assumptions, other than those made by Trivelpiece (1) 
have been made. For a detailed justification of these assumptions the 
reader is referred to Trivelpiece"s work. Extracting, the appropriate TE 
solution from Beck (9), the boundary value problem then becomes (referring 
to Fig. 1) 

Region I 0 ^ r ^ a (Electron beam within a plasma) 

From Appendix A , 
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E^AfVT'CT,*.) 



H z - A^T 0 (pyy) 



~ B& - f A%T, (p 

0 r 

Region II a £ r ^ b (Plasma region) 

From Appendix A 
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6 33 



E z = B,Hr a (T^ i BzWfcn.) 

t w p/i faTiJj(T*A.)tB i T 5 .V,(n.ft.jj 

Hi - A^To (p*) 

E s = (p'E) 

r 

Region II I b ^ r ^ c (Free space) 

Ez. = g J 0 (pj) + C^K 0 (f Jn -) 

c,i, (fiA.) - tm o %*-, (/> a.) 

r P 

Hz - c 3 2 o(p^) + fyKo (p' 1 ) 

B e = -t2j& C 3 I, fa) C y Kf ( P a) 



2.25 

) 



12 . 26 ) 



( 2 , 27 ) 



( 2 . 28 ) 
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Region IV r ^ c (Free space) 



Bz. ~ K 0 ( p a.) 

He =• -&^?D,K,(p^) 

P (2,29) 

H z - D$_K 0 (p».) 

P 

The boundary conditions at the helix are: 

ai -- at 

COT (2*30) 

Hi + Hi corf =Hi + HicoTY 

As the general problem is rather complex, several simpler cases will 
be considered first. These are 

Case I, the case of an electron beam completely filling the interior of the 
plasma and the helix, i.e. , a region I and IV problem. 

Case II, the case of an electron beam of radius a^, less than the diameter 
of the helix, with the plasma filling the helix, i.e, , a region I, II, and 
IV problem. 

Case III, the case of an electron beam of radius b ( passing through a plasma 
of radius b, surrounded by a free space region and a helix at radius c, 
with free space outside the helix, i.e. , a region I, III and IV problem. 

Case IV, the case of an electron beam of radius a_ passing through a plasma 
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of radius b surrounded by free space and a helix at radius c with free 
space outside the helix, i.e. , a region I, II, III, and IV problem 

Determinantal relationships are obtained for all four of these cases 
with the algebra relegated to Appendix B The determinantal relationships 
are given below. 

Case I 
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3. Method II, (Extension Of Brewer's Method) 

This method is an extension of the method of Rigrod and Lewis 6/ and 
Brewer (7). Attention is also invited to Beck's (9) excellent treatment from 
which this presentation proceeds. 

Let us consider an electron beam with no angular variation in charge 
density, fields or electron motion, i.e. , 




We may then write the Lorentz force equation in cylindrical coordinates as 

^ J ( 3 „ 1 ) 



[~v ^4 







4 



(3 3) 



provided we assume as does Beck and Brewer, that the contribution of 
E to the right side of (3.2) is small with respect to B r. More will be 

1 t7 Z 

said of this approximation later. Using Busch's theorem from electron 
optics, we may write 



m - 







(3.4) 



Where r is the initial radius of each electron and B~ the z-directed 
c 0 
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magnetic field at that point (assumed constant). This may be stated as 

0 — aT I W C 0 — I (3„5) 






Substituting this into (3.1), one obtains 



X = +f 



Let us now perturb (3.5) by letting 



_ 4 



e = & + 






(3,6) 



to obtain 






Jl~a 



<a>4 



'lo rc 






*0 






(3,7) 



Next, assuming that the DC ripple or scallop on the beam is small, let 



- -7^ + i 



Perturbing this expression as before yields 



^ 0 



8) 



ytj = ' r l ^ 



(3.9) 



Perturbing (3.2) and disregarding the effects of AC magnetic fields on 
electron motion compared with electric field effects 



1 



1 ( w+j'i'y 0 )^ 

Applying the continuity equation 



(3,10) 



17 






(3.11) 



v-x = 



to the AC charge density, one obtains 



P = A3 ViW 

' 0^p%V 




(3.12) 



At this point, let us examine what has been done. The equilibrium 
relation (3.8) has been perturbed to obtain (3 9). As stated in the intro- 
duction, electrons are assumed to interact only through the electric field 
and the above relations are derived as if the individual electrons compose 
a stream with continuity of charge in this stream maintained through 
equation (3.11). If a second group of charged particles were present, such 
that the DC equilibrium condition remained valid, comparable expressions 
for the time varying quantities could also be written for the second group 
of charged particles. The model of the plasma set forth in the introduction 
fits these conditions quite well since overall DC neutrality from the plasma 
ions and electrons is maintained. Assuming that the continuity equation 
holds (which it must as no mechanism for loss of charged particles has 
been assumed in the model), one can write expressions comparable to 

(3.9), (3.10), and (3.11) (with U =0) for both ions and electrons, but, 

o 

since the mass of the ions is so much greater than that of the electrons , 
their effects are extremely small at any frequency considerably above one 
megacycle per second and will be ignored. Designating the beam electrons 
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by subscript b and the plasma electrons by subscript a and dropping the 
subscript^ from field quantities since only ac fields are to be treated, the 
following relations are obtained. 






,ir (cv+j,Y(J 0 ) 






(3 o 13) 
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— r / ^ 

/it I n r 






(3 „ 14 









(3.15) 







(3.16) 
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, = v £ i- 
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(3.19) 
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(3 a 20} 
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^ Uu^-SiX 



Assuming that 




P 



XT 

oo, Iaxl 



-f- fi ajt ; 

/ n 



0-& /Sifr 



(3,21) 



and, from Maxwell's equations in a media containing charge 



t - X. Al± _ 4&/4 rr 



F 






(3.22) 



From (3.13), (3.14), (3.19) and (3.20) 



X = Pa.'f- Uj) lE«. -J- 

to*- J ?* ^ 



(3 .23) 



It will be convenient to write 



^ = fit 



A. 



(3.24) 



From (3.22) and (3.23) 

P — X£z. — F B 

pi 1 ^ (3,25) 

At this point, Beck shows that the second term of the last relationship is 
small with respect to the first for electron beams and may be ignored. 
This approximation is not as clearly well taken in this analysis and will 
not be made. Comment on the effect of making it will be made later. 
Solving for 
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Y is 



(3 26 . 



Small DC beam scalloping has already been assumed and since the 

assumption that the magnitude of the perturbation is small is inherent in 

the perturbational approach, r ^r, Further,, if we assume that the dc 

c 

magnetic field strength is everywhere constant and the cathode is not 

shielded, B v B and 
o z 
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* 



( 3 , 27 ) 



Making this approximation and writing 








(3 28 ) 



one obtains 





*%»• -l *'/$»?} 






( 3 , 29 ) 



Writing this as 




LE^~F 3 1 2. 




( 3,3 0 .? 



and 
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l£- 









(3 311 



and from Maxwell’s equations in a media containing charge 

1 1-U IL) -/>*£ = ^ 1 1LT) 

a dn, \ } a/ r *■ wO/6o ^6 0 a <?a( /y Y 



there results 



i i 



A ^A. ^ A 



^ $E Z - 0 



where 



v^ = _ P % (lf- M £o+L)(P' l U UI /'fi) 



32 ) 



33) 



34? 



The solution of this equation is, 



5 Z - ~l~C$ N 0 (y z si) 



(3 35) 



The notation at this point has become rather cumbersome. For clarity of 
presentation and ease of manipulation, some arbitrarily defined constants 
must be used. Although the use of the components of the dielectric tensors 
as derived in Appendix A, is objectionable in that an equivalent dielectric 
is not being used in this analysis, these quantities are familiar to workers 
in the field and their use will facilitate comparison of this method with the 
other two methods presented in this paper. Introducing these expressions 
and a quantity 
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(3 3 6) 



/f = A7 



the significance of which will appear later, we have 

f, = i w t<,(6r')+H(Prt,)e 0 



3?' 



L - -j,w£ 0 (£ 3 ~/ s ) 



(3.38) 




^ 6 . ( 6 ; ~ 6 „) 




3.39) 







(3.40) 



At this point, a convenient check on the development exists. If we 

let p . go to zero and make all the substitutions indicated in (3 34) 
ob 

we obtain 




(3.41) 



c 

The second term in the second bracket in the numerator can be shown to 
be absent if, in equation (3.22), the term containing J is ignored as 
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negligible. If this term is omitted, the result is then identical to that 

obtained by Method I (Trivelpiece's "slow wave r approximation), it is 

obvious that this term is not always negligible, even for very slow waves 

(ir’i >>k^). Similarly, if the term containing J is ignored in (3.22) 

and if we let p =0, the results can be shown to agree with Beck's 9 
oa 

confined beam development. 

Thus far, the TM solution has been obtained with the principal approxi- 
mation being that is negligible with respect to B^r-B^z in (3.2) 

This, in effect states that the coupling of the TE mode, from which E 

1 O 

is obtained, to the TM mode is negligible. An approximate TE solution 
which does not ignore the coupling of the TM to the TE mode will now be 
derived. From Maxwell's equations in a region containing charge we 
obtain 



A, 



+ mV)", 




(3.42) 




X. I T 



(3 43) 




and, from Appendix 1 



(3,44) 
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3 45 




CU^ — Utc*' 






-1 



-fj*(uj+j.x ^o) Eg. ’^su 

(Cu -t r Y^)^ - COE 



Assuming 



^OCL /V /&OL fo+ AS l&* r 



we obtain 






'^(AJufep 
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~yv 



00* — U ) ^ 

' ou c e 0 u& 

00 { * - uo a. 



(co+j.yu 0 y - lu£ 

co c £o 
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which, for brevity, will be written 

^9 = // (^i ~l) ^o. 



Substituting (3.44), (3.26), (3.36), into 3.49/, we ob am 





“/ 



A/4 



•f 



*fi*Wc.£o (tri) 3 5 , 
^ ^ 



We now have 



3 46 



{3 47 



3 48 



3 49 



3 50 
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— ft (£,-') 1 If JJ 






A. 



<5n. 






where 






(3.51) 



(3.52) 



Equivalently, let us write 

(3 o 53) 

where the primes denote differentiation with respect to r. 

A solution to this inhomogeneous Bessel equation may be obtained 
by letting 



r 



H " +f2HL - M.- - rM.t.a-0 

*- ~ ^ n. r z p.f 



// 






Ei + 



SL 



H z = A U,(^.) +■ 



(3.54) 




(3.55) 



Making these substitutions in (3.53) 
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(3.56) 



off -t- A -f U 

' SL 






-AK 



~h 



1A 



% 



Affr+CW&ieri) 



.3^ $ % 



r o 



For brevity, this will be written 

/(Vi) + {( L *») = 0 



(3.57) 



If we require that 




as we must to be consistent with (3.55), f(Uj) will then be identically 
zero. We may, without loss of generality, set B=C=1 . We then have 

— Z/f ^4 ~ 0 (3.59) 

the solution of which is 



Ux = ti 1 * (ft 71 ) +Ukj^) 



(3.60) 



which then furnishes the reason for the somewhat peculiar choice of 
constants made in (3.3 6) and 



A = F. = — 

* rM 



(3.61) 
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The complete solution is then 



+ (3 o 62) 

As was stated earlier, the plasma alone is considered as the limiting 

case of zero electron beam charge density. Hence we may write expressions 

for the plasma region alone which are identical in form to the preceding 
2 

with a.' = 0- The additional subscript jd will be added to distinguish 

between the two regions. 

As stated in the introduction, a distinguishing feature of this method 
of analysis is the replacement of the rippled beam with an equivalent 
cylindrical beam with a surface current density. Expressions for the 
surface current density will now be developed. Following Beck's procedure 
we write (at r=d) 





Using (3.15) and (3.16), we obtain 

^2 " U 0 ft " 6/) Ki (J) ‘ 3 • 65) 
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( 3 . 6 6 > 




€o 



MEJd) 



From Maxwell's equations, we obtain 




2<AJ6 0 

X 




+ _Lr 
T Y ^ 



(3.67) 










(3 . 68) 



Using (3.50), (3.24), (3.2 6) and manipulating the various constants 




(3.69) 




at, , 

t LM + 



( 6 t ~i) Q £~ 2 
p h 



(3,70) 





+ 




AM*) HMW 



(3.71) 
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(3 .72) 



£ e = + 




All the field quantities required for the boundary value problem have 
now been obtained. Now we shall establish the boundary conditions. 
From Stratton (11), we obtain the requirements that tangential E must be 
continuous and that 

/A X ~ //jr] ~~ ^ (3 o 73) 



which then provides the following set of boundary conditions 




The boundary value problem is then; 
Region I (beam and Plasma) 



(3 74) 




(3 o 75) 







(3.76) 
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H z - A { F^T 0 {X^) +A Z X 0 (fiA.) 



(3 77 ) 



e b = m w % fe-/) a, j, ai - m,x <n* 



n ~a 



+ 1 D 



% 



Region II (Plasma) 




Wz = b i^p^( t/ 1 ) +4 § P M> &/) + b 3 Jv ( p ¥ a) + # y /c 0 (p ¥ 





Region III (Free space) 

- C, 



(3 78 ) 



( 3 . 79 ) 
(3 = 80 ) 



) (3 = 81 ) 



(3 = 82 ) 
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( 3 , 83 ) 



(3 84} 



u - 40 th. 

'k /> 



Cil/(pyi) ~~C x K f (p^ 



H z - C- 6 I 0 (/)*.) + Q K 0 (/>^) 



£9 ~ ^ j/3-^/ *-</ (p' 1 ) 

Region IV (Free Space) 

£z - &i K (p*-) 

He = -*j k t>,K(p* 

H 2 - QsACH 

e© = ^5 K, (H 

As in Section two, the same four cases will be cons 



(3 . 85) 



(3.86) 



(3 „ 87) 



(3 88) 



(3.89) 



(3.90) 

with the details 



contained in Appendix C= It must be noted that, since this method calls 



for a rippled beam, space for the ripple must be allowed. Cases one and 



two must then be considered as limiting situations where the free space 
region has shrunk to zero. 
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4. Method III, Solution of Field Equations by Kales Method 



In Appendix A, tensor dielectric constants have been derived which 
take into account the effects of the electronic motion in an electron beam 
and plasma. Using the tensor dielectric constant, we may then solve 
Maxwell's equations 

V/ E ~ H (4, j) 



VX/J - T -t^coD 



V'B = 0 



(4 3) 



V'D - /O 



(4 4) 



for a region with zero charge and current density, since the effects of 
electronic motion are already taken into account by the tensor dielectric 
constant. It will be convenient to change the notation slightly so that 









0 


6 - 




e >< 


0 




0 


o 





0 



B-E 



(4 , 5) 
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the net effect of which is to include the factor £ 0 in each of the<^ 

Kales (8) developed a method for the solution of Maxwell's equations 
in anisotropic media such that the dielectric constant was isotropic but 
the permeability, ji, was a tensor of the same form as (4.5). This 
method has been used by Stafford (13) to study resonance phenomena in 
a plasma column, and by Johnson (14) for the case of a plasma in a cylin- 
drical drift tube. In the presentation that follows, Stafford's notation 
will be generally followed. 

Proceeding, we may write 



-h^ z ^(~Y)E 



(4.6) 




(4.7) 




(4.8) 




(4 



9) 



Separating transverse and longitudinal components 




(4.10) 
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(4 11) 



-VM =-*«K 

4x(-4 - wi) = r J6 "l. 

solving (4.11) for $ z / 

4*4 = + 4*5^ 

and noting that, if F =0 

z 

4Xf'S z X X) - Ov F& - (a z - a z )F =r - F 



we may obtain 



*> = - 4a w* 






(4 . 1 6) 



Substituting (4.14) and (4,15) into (4,13) and manipulating somewhat, 



we may obtain 

+ ^ VtE z - yg - ^ 



(4.17) 
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Taking the cross product of with both sides of (4.13) 

V t H z +Y — j. U) (q' l X 



and substituting (4.14) and (4.15) into (4,18), we obtain 







At this point, it will be convenient to define 

k * = 



(4.20) 



U/V 6 ^ 

(4 ,.21) 



from which follows 

Af CO 6)^— UJ = (f 



(4. 



Using (4.17) and (4.19) through (4,22), we obtain 

( K 1 / ,v ) ^ « 2 X tj. [- U‘\ t(^L. l)r f'%e„)£ z 



( 4 . 23 ] 









t K*m 
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Solving (4.11) for H , taking the cross product of a with both sides 

t z 

and substituting in (4.13), we obtain 

(4 24 

Substituting (4.11) into (4.18) and then using (4.24) and (4.20) through 
(4.22) we may obtain 



If the divergence of both sides of (4.25) is taken, the divergence of the 
second term on the right side of (4 25) can be shown to be zero by vector 
identities. We then have 



+ 0Jf4A'%) -f- 




(4,25) 



(K t A'' 1 ) V-B t = 



(4.26) 



Performing a similar operation on (4.23) yields 



yy) V-H k = \£( fwe^E z -K'Mi) 



(4, 27) 



We may write D = 6* E in the following manner 



D = 5, 



</ ^ Q Z €33 



(4.28) 
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If we take the divergence of D and note that (from (4.10)) 






we obtain 






o 





writing 



v-e = v-f*. - 



and making the necessary substitutions results in 



*•5- = - 




£m 

6 (/ 



Jr^. 



Let us now write 

V'H = o — v-ff t -y// z 



If we now substitute (4 . 32) and (4.33) in (4 . 26) and (4.27), we 
after some lengthy but routine algebra obtain 

N?X i-(K~A A p)H z + E z = 0 

^?£ z - r ^ E z ~ 0 

e n 



(4.29/ 



(4 30) 



(4 31) 



(4 32) 



(4 /33) 
may „ 



(4.34) 



(4.35) 
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For brevity, let these be written as 



V^Hz. ta.H z _ =0 



(4 3 6) 



V?E, + CE l tdfk = 0 



(4 37) 



In order to obtain a solution to this pair of simultaneous equations , 
the artiface of assuming that both and may be expressed as 

linear combinations of two other functions of r and 0 shall be used , 
The constraints between the constants that must exist to permit a solu- 
tion for the two new functions will then be determined. Having found 

these two functions, we may then solve for E and H and use the 

z z 

uniqueness theorem to state that these are the solutions. 

We shall now let 



Substituting these two relations in (4.36) and (4.37) and manipulating, 
we obtain 



£l +■ 



(4.38) 






(4.39) 



40 



(4.40) 



IT 1 ', + <4 



1% -r,l\ 



+ H}. 






= 0 



- %,( C A pjp ) 

Pfr- VP 



f W, 



/’/fef,r^)-(?.rc/;^,) 



= O 



(4.41) 



If we require the coefficient of in (4,40) and the coefficient of Uj 
in (4.41) to be zero, we will have two equations of the form 



tjV t -<? V = o 



(4.42) 



the solution of which may be written 



F = 



C,J m P>A) -PCjA/jAn) 



rh 



/ 



.0)9 



(4.43) 



Setting these coefficients equal to zero and manipulating we may obtain 



A 



* - c _n_ +*($1 



p, 



c !± tfl 

it 



(4 44) 
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2 - C/, 3 

^ n 



(4 45) 



Equations (4.38), (4,39), (4.44) and (4,45) comprise a system of six 
equations in eight unknowns leaving two relationships which we may 



specify arbitrarily , We could let p^=p^=l. However the choice 






a 






4 6) 



will result in a more compact notation. Making this choice, we find that 

A = a.+c±i(d- cf+tM ( 4 . 47 ) 

'A a 



where the upper sign is to be taken with subscript one. In terms of the 
dielectric tensor (4,5) 




It is also found that 




4(4 ~ c) 

ci 



(4,49) 



Before proceeding to find the field quantities themselves, it is to be 
noted that the method breaks down when the expression in the denominators 
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in (4.40) and (4.41) is zero. It can readily be shown that this occurs when 



(q-cf = -yjrj 



(4,50) 



which yields the following reguirement on ^ for this condition to exist 



y= 'fsl. ± 



~/r ^33 



4.51 



This condition will be discussed in Section five and for the present we 
shall only note that it occurs for certain particular combinations of the 
system parameters and that when it does we must examine the problem 
more closely. 

Using (4.38), (4.39), (4.46) and (4.49), we may write 

=£ 




4 .STMm tA IJVU.n)\£ 

4*0 *0 W A / 31 i ^ ' 












(4 „ 52) 




(4.53) 



as the most general expressions for E and H with the time and z 

z z 

dependence suppressed. We shall now limit our consideration to the 
axially symmetric case. 

Using (4.23) and (4.25), we may obtain the remaining field components. 
After some routine operations, we obtain 
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44 
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We may now proceed to solve the boundary value problem as was done 



in Appendix B for Method I. Since the procedure Is the same only the 
results will be given and the details dispensed with 

For case I, (beam and plasma filling a helix of radius c, we obtain 
equation (4.58) as the determinantal relationship Before proceding to the 



remaining cases it is convenient to define 

R = ^ 

-t l 



(4 59) 









(4.60) 



L-i 









/i 



.61. 



The additional subscript p will be attached to denote these quantities 
for the plasma region. Using these parameters, the determinantal relation- 
ships for cases II, III, and IV are given by equations (4.62), and (4.63) 
and (4.64), respectively. 
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5. Comparison of the Methods 



Trivelpiece's "slow wave" approximation is obviously the simplest of 
the three methods even though it is still so complicated that computer 
results are required for a clear insight into the system, a fact which de~ 
creases the value of its comparative simplicity greatly. It has several 
shortcomings. First, it neglects the ac magnetic fields in its basic 
assumption, which later has the effect of denying the possible existence 
of a TE type of solution. This leads to difficulties when attempting to 
match a field solution obtained from this method to any case where a TE 
solution is required, as at a helix. Equation (4.53) shows that a TE 
solution does exist and that it has poles, indicating that, regardless of 
how small the arbitrary constants, the fields associated with the TE mode 
will become appreciable in the vicinity of these poles; hence this approxi- 
mation cannot be used in these areas. One of these poles is a zero of 
as given by equation (A. 23) and indicates that the extension of Trivelpiece s 
method fails in an area not predicted by itself. Equation (3.41) also 
brings out a significant difference with method II. 

The extension of Brewer's method appears to improve the situation 
somewhat at a great increase in complexity. This method, by neglecting 
Eg, in equation (3.2), in effect, neglects the effect of the TE mode upon 
the TM mode but does predict a TE mode which is a function of the TM. 

The starting assumption in the analysis is that (neglecting the term 

in (3.2)) 
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S 1 



Using the relations developed within Section 2. this requires that the ratio 



be small. Examination of this relation shows that it has a second order 
pole at 



the values of Uj^ andd^ usually found in traveling wave tubes the ratio 
(5.2) will remain small. Also, the z term in 3 2) has been neglected 
since B^_ is very small. Looking at this term again shows that z has a 
pole at the frequency given by (5.3). The situation is no* clear and the 
determining factor is believed to be losses which have not been included 
in the analysis. Examining (5.2) further shows the possibilities for failure 
particularly for fast waves. However no simple statement of these condi- 
tions reveals itself. 

In general, this method is believed to be an improvement over Method I 
but, due primarily to the complex boundary conditions it leaves little if 
anything to recommend itself in preference to method III 

Method III, the Kales’ solution has several interesting features First 



5 ?_ ^f\ 

b^~ (r*+C) w(uj^ni) 




without corresponding zeroes. It is to be noted 
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coupled modes are predicted, i.e, ( TE and TM modes are related to each 



other by a constant and neither can be zero independently at other than 
special combinations of the system parameters,, Second, a two-fold mode 
degeneracy is predicted since we may take either of two values for the 
radial propagation constant* This aspect presents complications in the 
boundary value problem* 

Suppose that we try to match the boundary conditions using only one 
of the two possible modes. We would then find that only one arbitrary 
constant would appear in the field equations for region land matching 
boundary conditions would give us four equa tions in three unknowns . The 
problem would then be over specified and would allow solutions for, at 
most, particular combinations of the system parameters. By taking both 
modes of each degenerate pair, we may match any physically realizeable 
boundary conditions 

Very closely related to the above is the fact that the degenerate ) 
modes are not orthogonal. From Churchill (12), we may state the general 
requirements for orthogonality as 

a, U, (aS) -j- y{(cC) ~ o 
u, <*) + a^!u) -o isA) 

where U is a solution of a Sturm- Liouville system. Since these are 
precisely the relations which exist between and and H z and 

Eg for each of the modes individually, it is obvious that the azimuthal 
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modes (/V) ) are orthogonal for the same s., Since each of the s i is 

associated with a separate Sturm- Liouville system the degenerate (s.) 

1 

modes are not necessarily orthogonal. Evaluating an integral of the form 

Jq (/3-n.) dsi, (5 5 ) 

will demonstrate that, except for unusual combinations of the constants, 
the modes are not orthogonal. 

In summary, it may be stated that the attempted extension of a reia- 
tively simple method (Trivelpiece's) develops inconsistencies with the 
exact solution of the model and, until numerical results can prove to the 
contrary, must be used with caution. The extension of Brewer's method 
appears to improve the situation somewhat but the complexity of the 
boundary conditions deprives the method of any advantage of simplicity 
and hence, has little to recommend itself in favor of the Kales 1 solution 
The Kales' solution, while complicated in its derivation, is not more 
difficult to handle afterwards than any other. It has the distinct advantage 
of being an exact mathematical solution of the model, the only limiting 
assumption being the derivation of the dielectric tensor itself. It is 
readily demonstrated that the method fails at the poles and zeroes of the 
components of the dielectric tensor, and at a few other special combina- 
tions of the parameters such that the problem degenerates into a much 
simpler one. The failure at the singularities of the dielectric tensor is of 
no great improtance since the model also fails under these conditions, i.e. 
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infinite fields and propagation constants are not found in nature This 



method also has the advantage that it offers the possibility of extension 
to include the effect of collisions between charged particles since the 
conductivity of a plasma is a tensor quantity of the same form as the 
dielectric tensor. 
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APPENDIX A 



DERIVATION OF TENSOR DIELECTRIC CONSTANT 



It is permissible to treat a plasma as described in the introduction as 
an equivalent charge-free dielectric whose characteristics vary as a 
function of frequency. Further, the equivalent dielectric constant for a 
plasma in a magnetic field is a tensor quantity because the electric field 
vector and the displacement vector are no longer related by a single multi- 
plicative constant. The elements of this tensor are calculated by adding 
the convection current density to the free space displacement current 
density and setting the sum equal to the displacement current of the 
equivalent charge free region. Two cases shall be treated, first, that of 
a plasma alone, and second, that of a plasma with an electron beam 
through passing through it. 

PLASMA 

Neglecting a term as a second order quantity, one obtains 



Using the equation of motion and neglecting effects of AC magnetic 
fields on electron motion 




(A o 1 ) 




the following are then obtained. 
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-Z_ 

/») 










/YU) 



E /a -j-uf c AJ OL 



4 



Solving for the components of velocity yields 




^ = ~^ 

Substituting these components of velocity in equation (A.l) and 
for the tensor. 





6 // 




0 


6 - t5 




e n 


0 




0 


0 


^33 



(A„3 , 
(A., 4) 
(A - 5) 
(A, 6) 
(A, 7) 

(A. 8) 

solving 

(A „ 9) 
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Where 



£/, ~ I +* ouJl^ cu a . 




(A„ 3 0) 




J 



(A A 1) 





CU^ 



(A. 12) 



PLASMA AND ELECTRON BEAM 

Assuming that the interaction between beam electrons and plasma 
electrons takes place only through the electric field (in keeping with the 
model described in the introduction),, one may solve the force equation, 
and equation (A. 2) for the beam electrons and the plasma electrons separ- 
ately . Equations (A, 6), (A, 7), and (A. 8) are the solutions for the plasma 

electrons. For the beam electrons, assuming all ac quantities vary as 
jWt-ifz r 

e the force equation may be written 






(A., 13) 



from which, the following are obtained 




(A 14) 
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+ “'c'V 



Solving for the components of the ac velocity yields 



"k = ' 7 7 



(w^%Y-aJ^ ~f (WjtUoY-uj? 



(A .17) 



/ v* =_9 7 



-j.(un-$¥U 0 ) 


UJ C E,^ 


_(u 7+^j^Mjf 


(cu+jXy^-cu^ 



CA 0 I 8 ) 



Or — — "V) 



19) 



The convection current density is now written as the sum of the current 
densities due to the beam and plasma electrons and one then obtains 



+/&% =^g‘E, 



(A. 20) 



Solving for the components of the tensor 





6, 


J- 6 *- 


0 


of 

II 

Ut/I 


-* f 3. 




0 




0 


0 


% 



(A. 21) 
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where 



e - 1 , + 
'-I I ' 'a ^ — /iiSL / 



<4* 



w* -uj 9 ' 7 ^ c a -(ou+^ru , ) a 



6*) - 



U£ 

(A> 





w c 


1 

«* 

1 




pw^„) 


<V+#iW.)* 



6 ,= /- ^ 



UJ, 



m. 



uj^ ((MtyYUo)*' 



(A. 22} 



(A, 23) 



(A. 24) 
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APPEND IX B 



SOLUTION OF BOUNDARY VALUE PROBLEM , METHOD I 
CASE I (Region I and IV) 

Extracting the appropriate field components from Section two, we have 



Region I 




V A i VV> 


Wo !pr) 


^Kt^AjTjJjfTjr) 


£ -A, MU- 

2 p 


Region IV 




V B iV pr) 


H z =B 2 K 0 <Pr) 




E 

l P 1 


With the boundary conditions at the helix (r = 



^ = e£ 

< B - 3t 

A£ t /£cor 1 /' = Hzi-HgCory 

Matching these boundary conditions 

A^T 0 (r,c)~B,K o (pc) <b.4» 

~K ( f‘> o*) 
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A t »T t (T,c) = -coTf 



c)J 



^3. T o (/>*) + 



COTf 



*^-- 4 ,WV)] = 



B^K 0 (pc) 




(B- 7) 



Solving (B.4), (B,5), and (B.6) for A , B , and A 2 and substituting in 
(B.7); one obtains, after some manipulation 



J 

PCK>(PC) 



p3s 



K (ri 

KolPC) 



*>•7, (T,c) 



= 0 



(B„ 8) 



Making the approximation that 



and using equation (2,10), the determinantal relationship 



i | coT*-'rJt-i,(pc) 

pc K,(/>c) ~ fi- 



K, (pc) 

K 0 C/> c) 



507C) 



-0 (B 9) 



is obtained . 

CASE II (Region I, II, and IV) 

The appropriate field components, from section two, are 
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Region I 



H z - - 0 (p *^) 



"e = TfitTrf e 6 = -4 *f- 1, (H 



Region II 

Ez. “ 

** , *^[8,^37^ + ^W5*3 

Hz. — 0 ®^.) 

^5» ~ (p^l) 

Region IV 

^-'z — A 0°^) Wz ~ (p yx ) 

He =-[), if- K, (pd % = opf K) (r .) 



Ill 



(Bo 12) 



Matching tangential E and H components at the boundary between the 
beam and the plasma (r = a) 



CTf^) ~ B^^/i/ 0 (T^(L) 



(Bo 13 ) 



oja.) = [e,r a 3; (T^a) +b&h <*«)] 



14 



and eliminating the coefficients, one obtains 
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B, 



J o (J^a) _ £// 7a. 37 Ha Q) 



«/ 77^(7; a.) 



+R 



V.(%«i %uEM£h3i 



Jl ( TjCi ^ 6/Tj 1J (Tjci)~ 



- o (Bo 15) 



At the helix, boundary conditions (B,3) apply 

W r > c ) i-^wcnt) =o, /<„(/><) 




(Bo 1 7) 



£/ * X (Tic) + B^W t 0 » £ ) = 4 C0T ^ ^-X (P c ) 



A^IoCpc) tcor 




DxKJpc) -fycort^-K/fic) 



Solving the first three of these equations for , and D, and substituting 

in the fourth; one obtains, after some manipulation 
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B, 



Y^o ( 75. C )X 0 (pc) 1 Kq (pc) Y'TtfCJa.t) Cb7”V r ^ > '£ / / Tl'f'r r\ 

x, rr c) + KJrt ~Tp L WW 



- 



~h Bo 



'2 



^^o(Bc)T 0 (pc) /< 0 (pt) vy/j: A 
X,(pc) K, (pc) 1 0 * / 






- 0 



(B „ 2 0} 



Again making the approximation that 



one obtains 






-h 






tap 



(P C )K,(P C ) 



+ 



SSptffo#,, 6jJ w,^c)+ tiMmj 



(B.21) 

0 



Setting the determinant of coefficients of and B^ obtained from equations 
(B.15) and (B.21) equal to zero, the determinantal relationship is then 
obtained , 
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c 



-f 



-f jgtefzfcc) 



x 





K f « At fl»«) 




^(Vt) fe 


J~0 (1/ C P) 


v e< 63 o;(^) 




Jw r* 



6 / 37 (7Ta)_ 



-jJVqSML . CO tlA]. if, c s /\/.(t:c)+ KM) y n= £ >) r n 

r, v*>*, (p<) + a' ( ^ 0*) °" J ' '/J u 

CASE III (Region I, III, and IV) 

The appropriate field components, extracted from Section two 
Region I (Beam and Plasma) 



£* = 

«»= 



ttz. ~~ ( P ^ 

E e = *^1, (/>*■) 



Region II (Free space) 

= Ct^oCpn) iC)_/< 0 (p»-) 



n a = ^[gx,^-6^j] 

H-2 ~ ^3 (/W ~h Ko (p’ 1 -) 

Eq ~ ~ [^1/ (£ n ) ~ ^y Kj (p-^ 

Region IV (Free space) 

£z - £/ !<o (p^) Ht.- 

H a =-bpfzk l(H = ^KXpn) 



(B. 22) 



, are 



(B . 23) 



(B.24) 



(B. 25) 
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Matching tangential E and H components at the boundary between the 



beam and free space at r = b 

^X~C> + (f ko (pt) 



(Bo 2 6) 




(B 27' 



From the two equations immediately preceding, the following is obtained 





(P*) 
Io (/>*) 



{B, 28) 



h - ^3 



- Cu KAMI 

X, ( P 4) 



B 29 



It is then concluded that C. is zero and that A. = C„ . Matching E 

4 Z j z 

and H a 



4,n 0 ( t;*) = c,ijpt) +C&KJ/*) 



(B 30) 



* i yAU 7, (7-t) - f / X I M (p*) 



and eliminating , the following relation is obtained 



c, 



xm) r xpppj 

Xodi4) pt,T t q (JJ4) 



4 - 



K>M | ^K,(p4 ) 

X 0 (r,t) pm(U) 



-0 



(B. 32) 
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At the helix (r = c) » boundary conditions (B„3) apply 



C,X 0 Cpc) +C^ 0 (pc) =D,K 0 (pt) tB, 33) 

K/Cpt ) (B„34> 



C/X o Cp0 -hC x k 0 (pc) ~ Cj^Cpcjcorf 



(B.35 



C Z T 0 ( pc )tCoTr 






%Ko W + C0T * r-A ^ W ^ 



L 36) 



Eliminating D , D , and C„ , the following relation remains (after consi- 

1 L. O 

derable manipulation) 



C l 



(P c ) 


Jpcor'f 




A 

fc fpc) 


_ x > (pc) K, (pc) 


P a KtCpc) '} 


(/>£) 6 ^) 



o 



(B 37) 



Setting the determinant of coefficients of the C obtained from equations 

n 

(B.32) and (B„37) equal to zero, making the approximation that pZf-$ 
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and using the relations between the T and Y , the following determinants! 

n 

relation is obtained. 



Xo (pc) 


X.C/’*) 




I, (pc.) K, (pc-)_ 







Coc ) 

J i fa) K, (/><■) 



S0-*' 



-.♦^orV 

P r Ku? <■) 



Ko (pY) > l(i (p£) 

To 014) A /^3 



0 



(B ,38) 



CASE IV (Region I, II, III, and IV) 

The field components are as given in Section two for all regions. At 
r = a, the boundary between beam and plasma, the problem is the same as 
for Case II . The result is 





I 6 " ^ Ji£E2± 




1 e i <=3 Tift&J 



4“ 




At r = b, the boundary between plasma and free space 



&/ y c c 4~) ~ ~h k <>(/)&) 



(B.39) 



(B 40) 



Y 






(B.41) 



Substituting equation (B.39) in the above and defining 
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^ x 0 (r^mw qtwT.c Tjo.) 



'B 42 



H = A/ 0 t7^) J^a.) tytr^a )T, (r,«.) 



(B 43) 



After some manipulation, the following is obtained 



-L, 



(ft) -Jrh'sWH) -V*T, (ft](H(m)$ - 




>o <■ '0- 



'f £ 



X 



K 0 (p6) p\^/ t (T^b)G] -JpfyJHj 



■f 






At the helix (r = c) , the situation is the same as for Case III. We then 



have 



Tjpc) _ Ar^/^Oc) 

x, fy>c; />* *• 0»0 



X 6) 



K 

x, 00 



= 0 



(B 45) 



If the determinant of coefficients of the C obtained form equations (B.44) 

n 

and (Bo 45) is set equal to zero and the approximation that p^ 

is made and the definitions of the T used, the determinantal relationship 

n 

may be expressed as 
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/ 




APPENDIX C 



SOLUTION OF BOUNDARY VALUE PROBLEMS , METHOD II 
Case I (Region I and IV) 

Extracting the appropriate field quantities from Section three 
Region I (Beam and Plasma, 0 ^ r ^ c) 

(&c) 




*1 = AtF^TMA) *,(/j;c) 








Region IV (Free Space, r ^ c) 



f 2 - P, K 0 Cfc) 




- ^UOSrc 

P 



Dj /(, (pc) 



(ca) 



( 0 , 2 ) 



(C o 3) 



(C r 4 ) 



(C, 5) 



CO, 6) 
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Hz = c pz) 



(C 1) 




(C.8) 



The boundary conditions are derived from the basic conditions (3,74) 
and the boundary requirements of the "sheath helix" model. Brief comment 
has already been made that this case must be considered as the limiting 
situation of what amounts to case three (to be considered) with the free 
space region shrinking to zero thickness. If we view the problem in this 
light, the boundary conditions are found to be 




(C.9) 




(C ■ 1 0) 



E Zl = - COT r 



(C.ll) 



+ (%.i-^)corf = H Zjr 1-H 6jr CoTf 



(Co 12) 



Forming the quantities required by the last two relations, using (3,2 6), 
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(3.65), (3.66) and the relations given above for region I, we obtain 



- A I 



_ 






HsXo^c) 



(C. 13) 







+/?+* 




Applying the boundary conditions given above, we obtain 



D,KoC/=c) =A,T 0 OI h c) 



(C.3 4) 



(C. 15) 



(c .6) 



A K (f.0 = ^ COTtK, (pc) 

P 



(Co 1 7) 



4 / 

A 



W 



,-g 

p 



+A 3 J v^ c ) +■ 



-4 




/«{"$* l W / j 

Ak A> (/» 0 - A ^ (pQcorr 



COTf 



(Co 18) 



Using the same method as used in Appendix B, the determinantal relation- 
ship is found to be 
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/ //v ip fo-/) Jfc) tTCftO 






T 1 F) 



Ccrrf 



X o (ft, (■) + 



4M x > (K C J 

A 



*/?*■ 



(C. 19) 



-ft. +x , +A^(McoTr -h 

i w MPp\ M ) 



pZ(r*c)k,(fic) _j_ 

$<«/ta>T9 K, (pc) 



PKW 




It is believed obvious, at this point, that the complexity of this 
method is considerably greater than that of Method I. It is also believed 
that nothing is to be gained by writing lengthy, complex expressions from 
which no insight into the nature of the problem can be obtained by in- 
spection. The method of presentation of the determinantal relationships 
will, therefore, be changed to expressing the relationship in the form of a 
determinant set equal to zero. 

Case II. (Region I, II, and IV) 

Extracting the appropriate field quantities from Section three 
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Region I (Beam and Plasma, 0 ^ r ^ a) 

— A i ( )E/\/E) 






Afl -(/„ 




Hi ~ A, F x T c(^) i 'A 3 I l> (f /) ^) 



$ z t,(M -V f 

Region II (Plasma, a ^ r ^ c) 

Ez= B.TJtyi.) tB x K<~^) 







i w H 



I + 



i_ 

P* 

% 



WV 1 ) 



Hz. ~ 4 fifi Z ( «yi) V *6 ( / v a) FB^kJ n) 



(C, 20 ) 



(Co 21} 



(C . 22) 



(Co 23 ) 



(C„ 24 



(C. 25 ) 



(Co 26 ) 
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wcr 





(Co 27} 



Region IV (Free Space, r ^ c) 

fz = D>K 0 (/><■> 

K, (?<■) 

//z = ^7°^) 

5 * !±*jP K, (/>c) 



(C.28) 

(C. 29) 

(C o 3 0) 

(C.31) 



The boundary conditions for this case must be examined. At the 

helix (r = 0), the same considerations apply as were discussed in 

case I with the simplification that G is zero. At the boundary be- 

z 

tween the beam and the plasma, the situation is more complex. The 
boundary conditions expressed as equations (3.74) apply, but we must 
examine the expression for . We can write an expression simi- 

lar to (3.66) for the beam and also one for the rippled inner surface of 
the plasma region. Somewhat heuristically , the surface current 
arises due to the difference in the two media. Therefore, if we 
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let the difference between the two media approach zero in a smooth manner 
this surface current must also approach zero. must then be 




(C 32) 



After some manipulation, we obtain at r = a 






=4 




Hzrfy 






Xf 



+ 



(C.33) 



4 3 X„ (pa) + B, 7(1, a) + 

n *) ^ r / 

(C, 34) 

5 % 



and at r = c 



h z~ ~ 



^XOf^c) f- *&<*C€>(6 ir 0 Jf(Y C) 



Hi 



4~ (C.35) 



% 



"b d b 0 ( b/! -fi A/'ity C-) 



a# 



+ 



+fyK 0 (/%c) 
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By applying the appropriate boundary conditions at the interfaces 



between the regions, we may write eight equations in eight unknowns 
(the arbitrary constants) and write the determinantal equation as the con* 
dition for a non-trivial solution, i„e. , the determinant of the matrix of 
coefficients must equal zero. When this is done, equation (Co36) is 
obta ined . 

Case III (Region I, III, and IV) 

Region I (Beam and Plasma, 0 ^ r ^ a) 






(C.37) 



H e = i +j& 



, , r 1 



(C.33i 






(Co 39) 
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(C. 42) 



Hjp) = 4, J +$ i-r^l tfi) 

Region III (Free Space , a ^ r ^ c) 

£ z = Cl-^o Cp'O t* ^ Ko (P^) 








H z s £3 X& (/°*-) + K 0 Cp^) 




- ty k, (p^) 



Region IV (Free Space , r ^ c) 



E z = D, K„ Cp^) 






OUto 

p 




Hz. — &9- ^0 (/ 3a ) 



(Co 43) 



(Co 44) 



(C „ 45) 



(Co 46) 



(C» 47) 



(C . 48) 



(C , 49) 



80 



(C . 5 0} 



£ e = ^ D^H, (/>«■) 

Applying appropriate boundary conditions, as in Case II, the 
determinantal relationship is found to be equation (C.51). 

Case IV (Region I, II, III, and IV) 

Examining this case, we find that all the required quantities have 
been developed for the preceding cases, therefore equation (C.52) is given 
as the determinantal relationship without further comment. 
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